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We report measurements of the afo-plane penetration depth, A(T), in YBa2Cu307_{ films at various 
5. At optimal doping, critical fluctuation effects are absent, and 1/A 2 (T) from 4 K to 0.99 Tc is 
that of a clean, strong-coupling d-wave superconductor with Ao(0)/fcsTc — 3.3. As in crystals, 
underdoping reduces the superfluid density, n s (0) oc f/A 2 (0), without affecting the low-T slope of 
1/A 2 (T). These results, as well as electronic heat capacity data, are well described by an ad hoc 
model in which contributions to the superfluid and entropy are lost from regions of the Fermi surface 
occupied by the pseudogap. 
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A large body of experimental evidence indicates the 
opening of a k-dependent gap, or pseudogap, at a tem- 
perature, T* , above the superconducting transition tem- 
perature, Tc, in underdoped cuprates [Q. The pseudogap 
competes with the superconducting gap: T* and the frac- 
tion of the Fermi surface (FS) occupied by the pseudogap 
increase with underdoping, while Tc, the superfluid den- 
sity in the etfr-plane, 715(0), and the peak value of the elec- 
tronic specific heat coefficient, 7(T), at Tc decrease. A 
great deal of effort currently focuses on understanding the 
coexistence of these two gaps. Lee and coworkers propose 
that the fundamental physics lies in spin-charge separa- 
tion, a key element being the segmentation of the FS into 
regions either occupied or unoccupied by the pseudogap 
in the normal state [Q. With this in mind, we construct 
a simple model to describe our measurements of ns(T) 
and literature results for 7(T), in which only portions of 
the FS unoccupied by the pseudogap contribute to the 
superfluid and entropy in the superconducting state. 

We present new measurements of 1/A 2 (T) oc 715 (T) 
in YBa2Cu307_a (YBCO) films at various 6. We show 
that 1/ A 2 (T) in optimally- doped films is that of a clean, 
strong-coupling d-wave superconductor with a full FS. 
Underdoped films are well described by a Fermi liquid- 
like model, in which electronic properties are expressed 
as integrals over the FS, but the integrals extend only 
over sections of the FS not occupied by the pseudogap in 
the normal state. The fraction of the FS that survives 
is equal to the ratio of ns(0) of the underdoped film to 
ns t opt.(®) of the same film at optimal doping. 

The absence of critical fluctuations in optimally-doped 
YBCO films m and measurements of the effect of thermal 
phase fluctuations on 2D films of a conventional super- 
conductor lead us to conclude that fluctuation effects 
are weak in the underdoped films. However, the rela- 
tive importance of thermal phase fluctuations (TPF's) 
to single particle excitations is controversial. Carlson et 
al. II have shown numerically that a fluctuation driven 
superfluid density in Josephson junction ( JJ) arrays dis- 
plays features similar to some very clean YBCO crystals, 
namely, T-linear behavior at low-T M, Tc roughly pro- 



portional to ns(0), and a wide critical region |7Uq]. Tera- 
hertz measurements of the sheet conductance of BSCCO 
films also suggest a wide fluctuation region §. On the 
other hand, theoretical analyses of TPF's in underdoped 
cuprates conclude that fluctuations are too weak at low-T 
to account for the T-linear behavior of A (T) |[(J[nj. Con- 
sistent with these results, our estimates indicate that the 
effects of TPF's are minor. 

It is not known why fluctuation effects near Tc are 
weak in optimally-doped YBCO films || and some crys- 
tals 12 1, while appearing strong in other crystals MM. 



Films are of high quality, based on their T-linear A(T) 
at low-T and transition widths less than 1 K. Evidently, 
critical fluctuations are sensitive to structural differences 
which affect neither of these quality indicators. An esti- 
mate of the coupling between CuO bilayers in optimally- 
doped YBCO || finds that for T within 5 K of T c , the 
ratio of interlayer coupling to in-plane coupling (J'/J in 
ref. ||) is greater than unity and fluctuations should be 
strongly suppressed. To us, the presence of critical fluctu- 
ations in crystals is surprising, but their absence in films 
is not. We speculate that a significant portion of what 
appear to be critical fluctuations in crystals is, in fact, 
due to the rapid decrease in the quasiparticle scattering 
rate as T decreases below Tc [0 which serves to rapidly 
increase the conductivity, <j\(uj,T), for u) less than the 
gap frequency, Ao(T)/h, thereby increasing n$ || . The 
decrease in scattering rate is known to be less rapid in 
disordered samples jlq| . 

The d-wave theory which we use is an extension of the 
weak-coupling result for A 2 (0)/A 2 (T/T C0 ) JTtJ to strong 
coupling by increasing the ratio Ao(0)/fcsTco above its 
weak-coupling value of 2.14, while preserving the depen- 
dence of A o (T/T co )/A o (0) on T/T C0 - T co is the mean- 
field transition temperature. 1/A 2 (T) and j(T) are de- 
termined from the usual FS integrals over k-space and 
energy. 

For simplicity, we neglect possible deviations of A(k, T) 
from Ao(T)cos(20) which are reported in ARPES mea- 
surements on underdoped Bi2Sr2CaiCu20s+,5 (BSCCO) 
p8[ and which may or may not be present in YBCO. Pos- 
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sible anomalous behavior of the quasi- ID CuO chains @ 
is not included. Because such behavior is not observed 
in untwinned crystals, it is likely that chain specific ef- 
fects are negligible in our highly twinned films. Variation 
of the Fermi velocity, Vi?(k), with k is not included. If 
ns(0) is a factor F smaller than in the optimally-doped 
film, then FS integrals are taken only over the angular 
interval, ±F(7r/4), centered at each node in A(k, T). 

2D TPF's are predicted to suppress n s (T) @ by a 
factor: af Q (T)k B T^ X ± {T)/hR Q . 1/\±(T) = d/X 2 {T); 
d = 11.7 A for YBCO. The prefactor, a, is not well 
known theoretically. We use a = 0.175 which is consis- 
tent with measurements on a-MoGe films [|J and calcu- 
lations for a hexagonal array of resistively shunted JJ's 
|ll],[nj. R Q ee h/4e 2 ~ 1027 ft and f Q (0 < f Q < 1) 
represents quantum suppression of TPF's. The 2D tran- 
sition, T2D, is expected where X±(T2d)T2d — kTiRq /2k B 
= 9.8 mmK. Since T2D is near Tc this condition is 
approximately X±(T2d) — 9.8 mmK/Tc. (In the ab- 
sence of quantum effects, T2D is the Kosterlitz-Thouless- 
Berezinskii transition temperature Tktb-) As T ap- 
proaches T2D, the suppression of ns(T) grows rapidly 
due to nonlinear effects, reaching 20% to 50% just below 

T2D- 

Quantum suppression of TPF's in films has been con- 
troversial for some time. Recently, quantum effects were 
predicted jll| and observed ||] to suppress TPF's when 
ksT/h drops below the " R/L" frequency of the film. Ap- 
proximately, the sheet inductance is L = iiqX±(T) and 
the sheet conductance is 1/R — cti(w < Ao/ti,T)d. For 
optimally-doped YBCO, Jq should be much less than 
unity below 0.8 T c . 

Data presented here for films grown by coevaporation 
and sputtering are representative of other high quality 
films. Films allow the rapid and reversible adjustment 
of oxygen content without altering the thickness or mi- 
crostructure of the film. The coevaporated YBCO films 
were grown using the BaF2 method p2| with a room 
temperature SrTiOa substrate in an atmosphere of about 
5 x 10 -6 torr of O2, with a postanneal in wet oxygen. Af- 
ter careful refinement of the growth rates from measure- 
ments of film stoichiometry by Rutherford Backscatter- 
ing, YBCO films grown by this method consistently dis- 
play a linear low- T penetration depth. Additional films 
made by RF sputtering and by in-situ coevaporation of 
Sm, Ba, and Cu |^3| onto substrates held at 750 C to 
800 C are presented. Optimally-doped film transition 
widths were less than 1 K, based on the peak in the real 
conductivity <Ji(T). One of the codeposited YBCO films 
was deoxygenated three times to ~ 06.8, 06. 7, and 06.6 
(determined from Tc) using low temperature anneals in 
argon. 

The complex conductivity, a = g\ — ia2, is determined 
from the mutual inductance of coaxial coils driven at 50 
kHz located on opposite sides of the film [Q. With a 
well defined coil geometry and known film thickness, s, 
Maxwell's equations provide a as a function of the 
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The inverse of the 2D magnetic penetra- 



tion depth (1/Aj 



d/\ 2 oc n s ) for a YBCO film is 



shown with best low-T fits using a weak coupling gap 
Ao(0) = 2.14 ksTc (solid line) and strong coupling gap 
Ao(0) = 3.32 ksTc (dashed line). The data and the strong 
coupling fit are nearly indistinguishable. The upper right in- 
set shows four additional films with strong coupling mean-field 
fits A (0) = 3.2 ± 0.2k B T c (dashed lines). The lower left in- 
set compares the normalized electronic specific heat coefficient 
from Loram et al. [27] (open circles) with the d-wave result 
for A (0) = 3.32 k B T c (solid line). 

measured mutual inductance. Great care is taken to en- 
sure that measurements are made in the linear response 
regime to within 0.1 K of Tc by taking successive mea- 
surements at increasingly smaller drive coil currents. <J\ 
is very much smaller than 172 everywhere except very 
close to Tc- From cr 2 we define A^: 



A±(T) 



s/d 



(1) 



While A^(0) is determined to an accuracy of about 10%, 
limited by uncertainty in film thickness, relative changes 
induced by deoxygenation are known to better than 5%. 

Figure 1 shows l/Aj_(T) for an optimally-doped YBCO 
film. The solid curve is weak-coupling d-wave theory 
(A o (0)/k B T C0 = 2.14) fitted to the low-T data by ad- 
justing Tco- The dashed curve, which is nearly indis- 
tinguishable from the data, is strong-coupling d-wave 
theory with Ao(0)/fce = 3.3 TWo — 300 K, consistent 
with tunneling measurements |2q |. The upper right in- 
set to Figure 1 shows the excellent agreement between 
the measured 1/ Aj_ (T) and mean- field curves for several 
optimally-doped YBCO films and one film of 
SmBa2Cu307_5. Given 1/Aj_(0) and taking vp = 0.7 x 
10 5 m/s Hi, 7(T) calculated with the same gap com- 
pares well with data [^7) on bulk YBCO (lower left inset 
to Fig. 1, with 70 = 16 mJ/mole K 2 ). The small experi- 
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FIG. 2. The inset shows data for a coevaporated YBCO 
film at 3 stages of oxygen depletion, each successive stage 
having a lower 1/Aj_(0). Best low-T extrapolations are given 
by Ao(0) ~ 2.5 ksTc and are shown as dashed lines. The in- 
tersection between the data and the dotted line is where the 
classical 2D phase fluctuation transition would occur for com- 
pletely decoupled layers. The main figure displays the data 
from the first two deoxygenation steps and fits incorporating 
2D phase fluctuations (MF+TPF) which show a discontinu- 
ous drop at T2D- The upper solid line are the corresponding 
curves with fluctuation effects removed (TPF Corrected MF) . 

mental j(T) at low-T could be fit better if the model 
allowed a k-dependent vp and deviations of A(k) from 
cos(2</>). 

The inset to Figure 2 shows 1/X±(T) for a YBCO 
film at three stages of deoxygenation. It is striking that 
dns(T)/dT at low-T (not necessarily the normalized 
shape) is essentially independent of doping despite the 
introduction of oxygen vacancies into the CuO chains. In 
our model, the unchanged slope, dXj 1 (T) / cIT\t^q, indi- 
cates no change in the opening of the gap near the gap 
nodes [p8[ . The downward curvature of l/Aj_(T) reflects 
that of Ao(T) if TPF's arc negligible. Wc wish to extract 
the underlying mean-field behavior by estimating the ef- 
fect of 2D TPF's. For reference, the intersection of the 
dotted line with l/Aj_(T) (inset to Figure 2) would locate 
T2D were there no interlayer coupling and no quantum 
suppression of phase fluctuations. 

The dashed curves in the inset and main portion of 
Figure 2 are extrapolations of the data taken below T2D ■ 
The lower solid curves arc fits in which hypothetical mean- 
field behavior (upper solid curves) is suppressed by TPF's. 
Since the influence of TPF's is determined by T, the 
normal-state sheet resistance, and Aj_(T), all measured 
quantities, the only fitting parameter is Tco- The third 
deoxygenation step yields fits similar to the first two, but 
is left out for clarity. A reduction in fluctuations as a re- 
sult of interlayer coupling would bring the fiuctuation- 
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FIG. 3. n N (T) <x 1/A±(0) - 1/A_l(T) data for a YBCO 
film at optimal and 3 underdoped stages are shown. The data 
are truncated above T^d ■ To the right are the calculated nor- 
malized electronic specific heat coefficients for F=l, 0.8, 0.6, 
and 0.4 (from highest peak to lowest peak), k-space inte- 
grals in the model, are taken only over a fraction of the Fermi 
surface, F, as illustrated in the inset. 

corrected mean-field result closer to the data, so uncer- 
tainty in mean-field behavior is bracketed by the upper 
solid curves and the data. We take the simple extrap- 
olations (dashed curves) as reasonable approximations 
to mean- field behavior at least for T < T^d- Our con- 
clusions are insensitive to the extrapolation and to the 
"foot" seen just below Tc which is likely due to oxygen 
inhomogeneity. 

Figure 3 displays the normal fluid density, rijy(T) cx 
1/ Aj_ (0) — 1/ A j_ (T) , represented by the mean- field (dashed) 
curves of Fig. 2. Our model reproduces these when 
A o (0) ~ 300K is fixed and T co ^ T 2 / 5 , where F = 
ns(0)/ns. op t.{0) in the experimental range, 0.4 < F < 1. 
The gap ratio, Ao(0)//csT C 7o, increases with underdop- 
ing, but the gap ratio defined from the maximum gap on 
the contributing FS segments, A max — A (0)cos(f (1 — 
F)), increases only slightly. The contributing FS seg- 
ments and A max are shown pictorially in the inset to 
Figure 3. 

The model provides a basis for interpretation of j(T) 
pjjl . For optimal and mildly underdoped YBCO, there 
are three striking features in the data. One is that j(T) 
for T < 60 K is nearly independent of doping. Another is 
that the peak value of j(T) at Tc decreases much more 
rapidly than Tq with underdoping. Finally, the electronic 
entropy just above Tc, i.e., the integral of 7(T) from 
K to Tc, decreases significantly. This implies that the 
hypothetical normal-state 7(T) must decrease dramati- 
cally as T decreases below Tc, even though 7 is nearly 
constant for T > Tc- 

Figure 3 shows that j(T) calculated with parameters 
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fixed by A(T) has all of the experimental features men- 
tioned above. This is significant. Simple models which 
account for the reduction in ns(0) by an increase in the 
effective mass of electrons, for example, would not de- 
scribe 7 accurately. Our model incorrectly predicts that 
7(T) just above Tq should decrease with underdoping. 
We hypothesize that our model describes the electron de- 
grees of freedom that condense into the superconducting 
state, and that there exists, in addition, an anomalous 
contribution to j(T) which arises from degrees of free- 
dom not associated with the superfluid (perhaps from 
electron spin degrees of freedom) and which decreases 
rapidly in the vicinity of Tc- 

For strong underdoping, F < 0.4, the model predicts 
that the thermodynamic critical field, E>c{0), is propor- 
tional to nsiOy^Tco, and the upper critical field, Bc-M), 
perpendicular to the a6-plane is proportional to T co [g9| . 
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